% Ejemplo de solucion numerica de un sistema de ecuaciones 
% diferenciales ordinarias mediante el metodo de Runge-Kutta-Fehlberg.
% Requiere rkf45.m y arenstorf.m.

% Se ponen las condiciones iniciales
a = 0;
b = 500000;
ya = [12000;600];
M = 10000;    %h=(b-a)/M; 1000000
tol=1e-7;
clf;
hold on;
z= rk4('guanacos2',a,b,ya,M,tol); %Retorna [t, [x, y], h]
plot(z(2, :), z(3,:),"-;12000 guanacos, 600 pumas;");  %y(t)
ya = [12000;1000];
z= rk4('guanacos2',a,b,ya,M,tol); %Retorna [t, [x, y], h]
plot(z(2, :), z(3,:),"-;12000 guanacos, 1000 pumas;");  %y(t)
ya = [12000;2000];
z= rk4('guanacos2',a,b,ya,M,tol); %Retorna [t, [x, y], h]
plot(z(2, :), z(3,:),"-;12000 guanacos, 2000 pumas;");  %y(t)
ya = [12000;3000];
z= rk4('guanacos2',a,b,ya,M,tol); %Retorna [t, [x, y], h]
plot(z(2, :), z(3,:),"-;12000 guanacos, 3000 pumas;");  %y(t)
hold off;
xlabel('Guanacos')
ylabel('Pumas')


